Momentum classification of SU(n) spin chains using 
extended Young Tableaux 



Burkhard Scharfenberger^ and Martin Greiter^ 

^Institut fiir Theorie der Kondensierten Materie, KIT, Campus Siid, 
D-76128 Karlsruhe 

^Lehrstuhl fiir Theoretische Physik 1, Fakultat fiir Physik und Astronomie, 
Julius-Maximilians-Universitat, Am Hubland, D-97074 Wiirzburg 

Abstract. Obtaining eigenvalues of permutations acting on the product space of 
i-S^ ■ N representations of SU(n) usually involves either diagonalising their representation 

matrices on total-weight subspaces or decomposing their characters, which can be 
obtained from Frobenius' formula or via graphical methods using Young tableaux. For 
products of fundamental representations of SU(n), Schuricht and one of us proposed 
the method of extended Young Tableaux, which allows reading the eigenvalues of the 
cyclic permutation Cn directly off the, slightly modified, standard Young tableaux 
labelling an irreducible SU(n) representation. Here we generalise the method to all 
symmetric representations of SU(n), and show that Cjv eigenvalue computation based 
on extended Young tableaux is at least linearly faster than the standard methods 
mentioned. 
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1. Introduction 



Symmetries, wlietlier discrete or continuous, liave been a central concept in physics 
since its earliest days, and only by making use of them, explicitly or implicitly e.g. 
by choosing a suitable coordinate system, can most systems be treated analytically or 
even numerically. This is also true in the study of spin lattice models, a vibrant field 
of contemporary condensed matter physics that has produced many insights into novel 
states of matter and manifestations of order. A spin model consists of a cluster of 
spins arranged on some lattice tile, usually with periodic, but possibly other, boundary 
conditions, and a Hamiltonian describing the interaction of the spins with each other 
or with external fields. The 'spins' transform like some (irreducible) representation of 
SU(n), usually SU(2), and thus the Hamiltonian acts on a tensor product space whose 
dimension grows exponentially in N . In recent years, however, models with higher n 
have also received attention [Il[a[3ll2ElEl0IHll9l[TOl[IIl[Ta[ia[Tan5l[l6l[I71ttg 
especially since cold atoms in optical lattices hold the prospect of realising SU(n) models 
experimentally [HI [201 [2ll |22l [23] . 

The most commonly considered interaction is of the Heisenberg form H = 
J2i<j JijSiSj between spins on sites i and j with a coupling constant J^-. Such a 
Hamiltonian is inherently invariant under global SU(2) (SU(n)) rotations and conserves 
both Stot and Sl^^ (or, for general SU(n), highest total weight tUtot and total weight wl^^) 
respectively. Usually the Jjj obey some symmetry relations, often they even possess 
the full symmetry of the underlying lattice, implying that the Hamiltonian is conserved 
under all operations in =Sf , the (point) symmetry group of the lattice. Since in numerical 
studies the lattice is some finite tile containing sites, =Sf is a subgroup of 5*^?, the group 
of all permutation of N objects. 

In treating such a spin lattice model, either analytically or numerically, it seems 
clear that one should exploit the symmetries of the problem as far as possible. 
Therefore a product basis, where the z-components of all individual spins provides a 
complete labelling of all states, while conceptually simple, is not the best choice from a 
performance perspective. Rather, we should use a basis labelled by the eigenvalues of a 
maximal commuting subset of =Sf plus a number of other labels, e.g. the eigenvalues of 
as many further commuting permutations from Sn as are needed to provide a unique 
labelling. 

A mathematical problem that arises in this context is to determine the eigenvalues 
of these labelling lattice symmetries. Changing language from SU(2) to SU(n) the 
general problem can be stated like this: for arbitrary A^-fold product spaces V®^ of 
some irreducible representation (irrep) of SU(ra), find the eigenvalues the labelling 
symmetries L. The (n-l)-dimensional vector v is the highest weight of the irrep and 
has the same meaning for SU{n) as spin for SU(2). 

There are two traditional ways to solve this computationally, we know that the 
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tensor product decomposes into a direct sum of irreducible representations of SU(n): 

vr = © vT- (1) 

w 

Say we want to know the eigenvalues for our labelling lattice symmetries on the 
subspace V^"'" of all irreps Ko of SU(n) {i.e., appears a.^, times in V^®^). Then 
one way is to use character theory to determine the irreducible representations of 
Sj\f contained in V^'^'^. The eigenvalues of the permutations are then obtained by 
diagonalising their representation matrices in these irreps, which is possible for all of 
them simultaneously since they commute. Alternatively, we could simply write down 
all product states (p^ ^ Kf^' which have the total weight w^^^ = w (in the case of 
SU(2) this corresponds to a fixed S^^^ subspace), determine the representation matrices 
of the labelling symmetries and again diagonalise them (all simultaneously). In both 
cases, we need to repeat the process for those highest weight multiplets w' > w, that 
are contained in the subspace of the total weight w^^^ = w. In the case of SU(2) for 
instance, it is sufficient to consider the next higher 5*^^^ subspace, i.e., S^^^ = S+1, and 
disregard all sets of eigenvalues which appear in both subspaces. 

Young tableaux are a diagrammatic technique originally invented to compute 
various properties of irreps of the permutation group Sn [24J, but they and other 
techniques based on them have since seen a myriad of uses in both mathematics [251 126] 
and physics [271 EH]. Their usefulness is mostly related to the Schur-Weyl duality, which 
makes a direct connection between irreps of GL(n) and Sn contained in tensor products 
of some elementary GL(n) irrep. Due to this, YT also provide an elegant means of 
obtaining the decomposition of V.^'^ into irreducible representations of SU(?t,), i.e., of 
obtaining the in (JT]). 

A further use of YT was introduced in a 2007 paper [29] (see also [30j), by Schuricht 
and one of us: the method of extended Young tableaux. It allows the spinon content of 
an eigenstate of the Haldane-Shastry-Model for a chain of fundamental SU(n) spins to 
be read off directly from slightly modified Young tableaux. The interesting consequence 
we want to point out here, is that this also allows one to find the state labels \wtot,p) 
for the symmetry group of the ID chain. This group is the cyclic group generated 
by the single-site translation Cn where each eigenvalue 7 of Cat can be identified with a 
momentum p along the chain via the equation 7 = exp[ip]. To our knowledge, this is the 
only such method working directly with Young tableaux and as a practical consequence, 
this enables a significant speed-up of CTv-eigenvalue computations. 

In this paper we generalise the method of extended Young tableaux to higher 
representations of SU(n) and present numerical evidence that it indeed gives the correct 
results for the eigenvalues of Cat, as long as the representations which are coupled are 
symmetric, i. e., correspond to Young tableaux with a single row. We also find, that while 
it does give the correct distribution of momenta on the subspace V^""^ of all multiplets 
A, it does not assign these momenta to the individual YT in a way that would allow 
deducing the irreducible S'Tv-representation content of V®^^. A positive result would 
both have been quite useful in itself and would also have given a physical meaning to 
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individual YT, similar in spirit to the connection between YT and angular momentum 
states sought by McAven and Schlesinger |3T]. Lastly we show how working directly 
with the Young tableaux of an irreducible representation speeds up the computation of 
Cat eigenvalues over traditional methods by at least a linear factor A^. 

This paper is organized as follows: In section 2 we briefly review the method 
of Young tableaux for and how it relates to the irreducible representations of 
SU(?7,). In section 3, we restate the extension procedure for fundamental representations 
and reformulate it in a way better suited to both numerical implementation and 
generalisation to product spaces of higher representations. We complete this task in 
section 4, which contains the main result, namely that the procedure appears to work 
for all symmetric higher SU(n) representations. Furthermore, we comment on relation 
between the extended Young tableaux method and irreducible representations of Sn 
contained in the product spaces. In section 5 we compare the computational complexity 
of the extension procedure to that of generic methods of obtaining the eigenvalues of 
Cn- Finally, we summarize our results in section 6. 

2. Young Tableaux and SU(n) 

We begin with a short summary of Young tableaux and some of their traditional uses. 
A Young diagram or shape is a graphical depiction of an integer partition 



as k left-justified rows of boxes, where row j has Xj boxes in it. A Young tableau (YT) 
on the shape A is any filling of the boxes with integers between 1 . . . A^ (see Fig. [TK and 
lb ). Counting Young tableaux subject to certain building rule is what is at the heart 
of their application in the representation theory of both the symmetric group S^, the 
group of all permutations of A^ distinguishable things, and SU(n), the group of special 
unimodular, complex n x n matrices. 

Sn- — The irreducible representations of the symmetric group Sn can be labelled 
by integer partitions A, |A| = A^. Originally, Young tableaux where invented to provide 
a graphical method of computing the character x'^(-P) of an element P G Sn in the 
irreducible representation A [2^. Of special significance is the character of the identity, 
since it is equal the dimension of an irreducible representation: x'^(id) = dim(A) =: A^a- 
This dimension can be determined by counting all standard Young tableaux on the shape 
A. A tableau is called standard if the numbers in its boxes are strictly increasing in both 
rows and columns. FigHfc) for instance shows all standard YT on A = (4, 2). There are 
9 of them, therefore the representation (4, 2) of Sq must be 9 dimensional. 

Special unitary group SU(n) — In the context of SU(n) Young tableaux can 
be applied to decompose tensor products of representations. An irrep of SU(n) is 
characterized by its highest weight, where a weight is the (?T,-l)-dimensional vector of 



(A) = (Ai, As, ... , Afc), Ai > A2 > . . . > Afc > 0, 
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Figure 1. (a) the Young diagram or stiape to the partition (5,3,2,1), (b) the same 
diagram as a (general) Young tableaux (c) all standard Young tableaux on the shape 

(4,2) 



eigenvalues of the simultaneously diagonalisable group generators spanning the [n-l)- 
dimensional Cartan sub-algebra of su(n), the generating Lie algebra of SU(?2). A weight 
w = {wi, W2, . . . , Wn-i) is said to be higher than a weight w' = {w[, . . . , if the 

first nonzero entry in w — w' is positive. In the well known case of SU(2), for instance, 
irreducible representations are characterized by their spin S, which can be integer or 
half- integer (S'=|,l,|,...), and the highest weight corresponds simply to the highest 
possible value of S^, which is = S. Since SU(n) is defined as a matrix group, 
one representation is always the group itself, its carrier vector space being C". It is 
irreducible and of special interest, because by forming tensor products of multiple F„ 
and projecting onto subspaces of appropriate symmetry we can form all irreducible 
representations of SU(r;,), This is the main consequence of the the already mentioned 
Schur-Weyl-duality. 

Thus, Schur-Weyl duality in effect implies that one can use Young tableaux to 
decompose tensor products of F„ (or indeed higher SU(?t,) representations). Associating 
F„ with a single box Young diagram, there is a neat diagrammatic way to do this: we 
simply construct all A^-box standard Young tableaux with no more than n rows, which 
can best be done using the branching rule (see e.g. [25] )• The process is illustrated in 
Figure m This also means there is a 1-1 correspondance between an n — 1 dimensinoal 
heighest weight w and an n— 1-row shape A, so we can from now on use shapes not only 
to index S^- but also SU(r;,)-irreps. 

It is possible to generalise this procedure to decompose the product space V^®^ of 
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Figure 2. Building higher irreducible representations from the fundamental one via 
branching in the case of SU(2) and SU(3). Counting dimensions as a check, we see that 
the decompositions are complete in both cases: 2'' = 2 + 2 + 4 and 3'^ = 1 + 8 + 8 + 10. 
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Figure 3. Decomposition of a 3 x (5 = 1) product space of SU(2) 

arbitrary irreducible representations associated with the shape a = (cti, . . . ,ai). In the 
case of single-row a {^symmetric representation\ I = 1), which is all we will need in this 
paper, this generalisation is straightforward: We again use the branching rule, i.e., add 
boxes step by step, but now each number j = 1, . . . ,N appears \a\ times instead of only 
once and we have to take care not to put two boxes with the same number on top of 
each other (see Fig. [3]). 

As we have stated, the number Nx of standard YT on the shape A equals both the 
dimension of the irreducible representation labelled with A and the multiplicity of 
the irreducible SU(n) representation associated with A (via the correspondence highest 
weight ^ Young diagram) tensor product -F®^ of the fundamental representation of 
SU(r;,). This just another consequence of the Schur-Weyl-duality: one can show that 
the subspace V^^ in the tensor product F^"' always forms an irreducible representation 
of the permutation group Sn equivalent to the irrep labelled by the integer partition A. 

For an example we can again consider Figure |2l the two YT of shape (2, 1) tell us, 
as mentioned, that F^^ contains two doublets but also, that the states of these doublets 
(for each fixed S^^^) transform like a standard representation (2,1) under the action of 
the group 5*3. 

3. Extended Young tableaux 

We now turn to the problem described in the introduction of computing the eigenvalues 



of the cyclic permutation C^v on total spin subspaces K 



Let us first review the 



extension procedure for YT of fundamental representations introduced in [29] . 

Rule, original version — Let T be a standard YT of size A^. By sliding them to the 
right where necessary, arrange all boxes of T such that in each column of the resulting 
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Figure 4. A simple example of the extension procedure: the lower row slides to the 
right, s.t. '2' is above '3'. How many dots are placed depends on the SU(n) under 
consideration 
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Figure 5. Building extended tableaux box-by-box. Since we consider SU(3), three 
rows are marked with dots. However, the momenta assigned do not depend on n 

extended tableau the numbers in the boxes are in sequence {i.e., i above i+1 above i+2 
etc.). This will often require leaving empty spaces between boxes (see Fig. H]). Mark 
each by a dot. To each dot i we assign a number Oj in such a way that the average of 
all ttj within one column equals the average of the numbers in all boxes in that column, 
where the have integer or half-integer values with a spacing of 1 between the numbers 
from one column. 

HaldaneShastry model. — This version of the rule betrays the origin of the 
procedure: it comes from the physical problem of the Haldane-Shastry spin chain [321 
|33] . which consists of spins on a circle with a Heisenberg-type JijSiSj interaction 
where the coupling Jij = |?7i — ?7j|~^ decreases quadratically in the chord distance. In 
the original model the Si where S = 1/2 SU(2) spins, but it has been generalised to 
fundamental irreps of SU(n) |H [5l IMl EB]. This fully integrable model [36] has a singlet 
ground state and the excitations are spinous, which can be thought of as delocalised 
domain walls ('half a spin flip') in a background liquid with strong antiferromagnetic 
short range correlations. By interpreting each dot in the extended Young tableaux 
as a spinon, the eYT allow obtaining the spinon content of the eigenstates (which also 
conserve total spin and total momentum) and moreover assign each spinon a momentum 
number pi connected to the from above via 

Pi = 271 (ai - 1/2) /N. (3) 

The total momentum of a state is obtained by summing over all individual spinon 
momenta while the energy is essentially the sum of the squares [29]. In both quantities 
we need to include the constant offset po and p^ respectively given by 

7X — 1 

Po = Tx N (4) 

n 

The original rule is a succinct formulation of the basic idea, but it is not well suited 
to implementation on a computer, nor does it generalise directly to higher SU(?2) irreps. 
If one wishes to use extended YT for computing eigenvalues of C^r on tensor product 
spaces of higher SU(r;,) irreps, it is better to make use of the branching rule. Let us 
therefore reformulate the building rule. 

Rule, new version — Given a standard YT T, we start with an incomplete extended 
tableaux Ei{T) containing only the single box labelled '1'. we then add a box labelled 
'2' to El by looking whether in T, '2' appears in the first or second row. In the former 
case we add '2' to the right of '1', while in the latter case we put it below '1'. This 
yields an, i.g., still incomplete, 2-box extended tableaux which we call E2. We go on 
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Figure 6. The complete list of extended Young tableaux for TV = 4 x (S' = 1/2) spins 
of SU(2), rriT is the integer momentum number Apx/'iTT. 



building the full extended tableau E(T) step-by-step. In step k, having constructed the 
extended tableaux Ek, we obtain the next one, Ek+i, in a similar way as E2: we look up 
in T the row index r^+i of the box 'k-|-l' and compare it to the one of 'k', which is r^. 
If now Vk+i < Tfc, we add a new column at the right side of Ek and put the box 'k+1' 
in its r^+ith row. Otherwise, i.e., if r^+i > r^, we add 'k+1' into same column as 'k', 
also in the r^+ith row. The resulting extended tableaux we call Ek+i{T). After a total 
of — 1 additions we thus arrive at the final tableaux E{T) = En{T). 

This procedure is now directly implementable and above all, generalises to products 
of higher (symmetric) representations. To compute the total momentum pt associated 
to the tableaux T, we combine the spinon momentum numbers Oj in each column c of 
the extended tableau E{T) into one column number be 

= E - 1) (5) 

Written directly in terms of the average of the box-labels (i)^ and the number of boxes 
kc in the column, the be are: 

be = {n-ke) 1/2) (6) 

the momentum pt associated with T is then simply 

The momentum offset 60 •= ~ l)A^^/2 is still necessary to ensure that the sum ([7]) 
is a multiple of n. Thus no matter which SU(n) a tableau T pertains to (although 
clearly one must haven > 7^ rows in T), it is always assigned the same momentum by our 
procedure. 
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Figure 7. Higher representation extended YT can be built box-wise as well. 
Ambiguity below which of two identical fc-boxes to place a (A; + l)-box requires an 
additional rule 



In Figj6]we show as an example the extension procedure for all total spin multiplets 
of the tensor product space (1)^"^ of four S = 1/2. Since this is the tensor product of a 
fundamental representation, the shapes of the YT immediately tell us what irreps of 5*4 
the total spin multiplets belong to. It is thus easily verified, that our procedure gives 
the correct values: The lone quintet 5**°* = 2 must be fully symmetric and has therefore 
momentum 0, the three triplets 5**°* = 1 form the standard representation of 6*4 (of 
dimension 3 and associated with the partition (3, 1)) while the two singlets belong to 
the self- conjugate irrep (2,2) (2 dimensional). 

The Haldane-Shastry model is valid for Si G SU{n) not only for n = 2 and the 
mechanism of constructing excitations remains the same. Therefore this connection of 
extended YT and HSM eigenstates exists not only for SU(2) but higher n as well and 
since extended YT correctly describe the eigenstates of the Haldane- Shastry model in 
all cases, it provides the strongest argument in favour of the correctness of our procedure. 
A rigorous mathematical proof would still be desirable however. 

4. Higher Representations 

While the spins in the Haldane-Shastry model can transform like any fundamental 
SU(?T,) representations, not just SU(2), a straightforward generalisation to higher 
representations runs into difficulties. One can nonetheless hope, that the mathematical 
statement remains valid for a suitable generalisation of the rule, which is what we 
investigated. It turns out, that for products of symmetric SU(n) representations 
(represented by single row YT with \ {a) \ = a boxes) the box-by-box approach to building 
extended tableaux generalises almost directly (see Fig. [Tj). 

We have to introduce only one additional condition coming from an ambiguity in 
where to put a box ^k+V if there are several eligible open columns with boxes 'fc'. Given 
a tableau T representing a multiplet in {i.e., each number appears \a\ times) we 
require that the resulting extended tableau E(T) is minimal, i.e., has as few dots/empty 
spaces as possible. The way to achieve this is to consider the boxes with number ^k+V 
in increasing order of their row index in T, i.e., place higher boxes first. 

Consider the fourth YT from the top in Fig. [8l Assume we have already placed all 
'I's and '2's, giving us an (incomplete) extended tableaux with three columns, the first 
containing '1' above '2', in the second and third a '1' and a '2' respectively in the first 
row. If we now were to place the '3' from the lowest row first, then we have two eligible 
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columns (the first and third). Depending on our decision, we would end up with two 
different extended tableaux, one with momentum p = the other with p = 1. Placing 
the higher '2' (the one from the second row) first, there is no ambiguity, and we identify 
the extended YT with p = 1 as the correct minimal one. 

The only further change is a that the momentum offset bo acquires a factor |cr|: 

bo{a) = -^N\a\ (8) 

The \a\ dependence arises naturally as is explained in Appendix C. 

Figures [H] and M show examples of the extension in the case of |cr| = 2. In the 
absence of a rigorous mathematical understanding why it works and without the physical 
interpretation that backs up the extended YT procedure in the case of fundamental 
SU(r;,) representations, we checked the statement numerically for several SU(n) up to 
n = 4 and = 16 and find it does give the correct eigenvalue distributions (see Tabled]). 

Figure M illustrates a hmitation of our method. As mentioned, multiplet subspaces 
^®a,A comprise in general more than one irreducible representation of S]y. 

Only for products of fundamental representations {\a\ = 1) do irreps of Sn and SU(n) 
coincide. The as=2 = 6-dimensional subspace of all quintets in {S = 1)®*^ for instance 
{i.e., 1^1=^) contains the three irreps (3,1), (2,2) and (4). If we try an ad-hoc 
identification of these irreps with the YTs based on the latter's structure (see Fig. |9]), 
we see that the momenta assigned by our procedure and the momenta one would expect 
from this identification do not match: the three tableaux we identified as belonging to 
irrep (3, 1) are assigned momenta 2,3 and 0, while one would expect 1,2 and 3. 

Thus, while extended YT do produce the correct frequencies of momenta for each 
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Figure 8. Here are some of the extended YT we find when combining = 3 rep 6 
(=(2, 0)) of SU(3). The last column gives the tally of momenta for the shape. 
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Figure 9. The complete decomposition oi V = S = 1 for ^tot > 2 and momenta 
assigned by our method. The right column shows the ad-hoc identification of YT with 
the irreps of S4 as described in the text. 



subspace y^"^ as a whole, they give no help in identifying the irrep content 
of multiplet subspaces (beyond what the momentum frequencies themselves already 
reveal) . 

5. Fast tableaux generation 

In this section we want to elaborate how working directly with the Young tableaux, 
allows a useful speed-up of C^v eigenvalue computations in the limit of large and n. 
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Table 1. Unitary groups SU(n) and maximal tensor powers N for which we verified 
the correctness of extended Young tableaux 
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In the introduction, we already mentioned briefly two traditional computational 
methods for obtaining the eigenvalues of the cyclic subgroup generator C^. They are 
character theory and diagonalisation of the matrix of Cat on total weight representations. 
Both begin by writing down a product-state basis of a total weight subspace 
ly^ot = w. In the case of {S = 1)®'^ of SU(2) for instance this would be the 10 dimensional 
space spanned by 

Q|l,l,l,-l),Ci|l,l,0,0),Cl|l,0,LO) (9) 

where i,j = 0..3, k = 0,1 and the cyclic permutation C4 is applied to a state in the 
natural way. Clearly, these states form the basis of a representation of S4. 

Character theory -There are in fact two methods based on group characters: one 
working with the characters of Sn and another, simpler one, using the characters of 

The former, mentioned here for the sake of completeness, obtains the Sn character 
Xw of the total weight representation spanned by (referred to it simply as the 
representation from now on), i.e., we compute the trace of the representation matrix 
of one element from each conjugacy class in Sn and then decompose this compound 
character using the formula ax = '^[p] \ \.-^]\Xx{P)Xw{P) , where the sum runs over all 
conjugacy classes [P] C Sn where P is some representative of the class and \ [P]\ is its 
cardinality. The eigenvalues of Cjv follow directly, since each irrep A comes with a fixed 
set of eigenvalues. 

One arrives at more efficient way of using characters by realising that is also 
a representation of which means we can apply character decomposition directly to 
the representation matrices of Cn, Cjf, . . . , C^. We can thus compute the multiplicity 
fm of a momentum number m via the group characters of ^jv: 

TV 



■^"^ ~ AT S 



k=l 



271 i , 
mk 

N 



TrC^ (10) 



where m = 0, . . . , N — 1 labels the irreducible representations and k = 1, . . . ,N the 
classes in and Tr is the trace of the x n^, representation matrix of the kth 
power of Cjv- 

This is both faster than the full character decomposition and does not assume 
prior knowledge of all the irreducible Sn characters (which would in practice have to be 
computed too). We do however have to generate all the powers of Cn, which takes 
(at least) 0{N Uw) steps. 

We are not done yet however, for remember that the total weight representation 
iSS^ contains not only the irrep w but also some with highest weight w' > w, which we 
need to sift out. The SU(2) case is straightforward: we simply run the procedure twice, 
once for S^^^ — S and once for S^^^ — S-\-l and then subtract the -S'iv-inomentum tally of 
the latter from that of the former. The general case requires more work however: first, 
we need to know the positive integers c^w' recording how many states a representation 
w' contributes to w. The fact that only w' > w contribute and w is contained exactly 
once means that viewed as a matrix, (c) will be upper triangular with only Is on the 
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Figure 10. (Colour online) The standard YT branching graph of the shape A = 
(2, 2, 1). For all depth j the arrow-labels denote into which line the index N—j + 1 is 
to be put. Thus, the paths through the BG correspond 1-1 to all standard YT. An 
example of such a path and the YT it corresponds to is shown in red. 



main diagonal. To obtain the tally of all momenta for 11;*°*, we then need to take linear 
combinations of a certain number of rows of this matrix, such that, all contributions 
of higher multiplets are cancelled. 

Thus, total asymptotic complexity is 



This is however still not as good as the conceptually simple diagonalisation we will turn 
to next. 

Diagonalisation. — Diagonalisation is straightforward: we write out the representa- 
tion matrix of Cn in the basis and diagonalise it. In general, diagonalisation is of 
(time) complexity O {m^) for an m x m matirx, but since we are dealing with permu- 
tation matrices (in each row and column all entries are zero except for exactly one '1'), 
O (m) steps suffice. Like with the previous character methods, we will also obtain Cn 
eigenvalues belonging to SU(n) irreps of higher highest weight which can be got rid of 
in the same way, incurring the same factor. 

In the end therefore, diagonalisation is faster than the based character methods 
and if we assume that the representation matrix of Cn can be written down in O^n^) 
steps the final time and space requirements are 



Extension. — The extension procedure on the other hand works directly with the Nx 
Young tableaux on a shape A, assigning each a momentum number m{T) = 0, . . . , A^ — 1. 
The key to making it superior to the other methods, is that it is possible to combine 
YT creation, extension and momentum computation efficiently into one procedure. 

Let us first consider YT generation: one should exploit the branching property 
somehow, but a naive ansatz building up tableaux by adding box after box starting from 
scratch for each tableaux will require O {N Nx), which will only be marginally better 
than diagonalisation in the most interesting cases (the states with low total highest 
weight, e.g. SU(n) singlets) and due to the more intricate nature of the algorithms 
involved probably turn out to be somewhat slower in the less interesting ones (states 
with total highest weight close to the completely symmetric one). 






(12) 
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We can achieve a complexity of O (A^a) however, if we store the branching 
information in a suitable way: the branching graph (BG). It encodes the relations 
between a shape A and all (valid) shapes /i C A derivable from it by repeated regular 
removal of elementary shapes a in the form of a directed graph, where a was the shape 
associated with the SU(?7,) irrep from which we build our product space. 

The nodes of the graph are the shapes yU C A (with A being the root) and a labelled 
edge (/i —7- v\ I) goes from shape /i to if and only if the latter can be obtained from 
the former by a regular removal of one elementary shape a. A regular removal is the 
inverse of a regular addition, which is defined as the addition of |(j| boxes such that the 
resulting tableau is valid. The label / will be a list of length \a\, recording into which 
row we put the first, second, third,. . . , |cr|th box. An example of a branching graph 
for the standard YT on shape (2,2, 1) is depicted in Fig. [TOl Since |cr| = 1, the label 
consists of a single row-index only. 

As long as (J is a single-row tableau, as we always assume here, there will be at most 
one edge between nodes. However, the branching graph is also defined for tableaux built 
from multi-row-cr (non-symmetric SU(n) representations), but there it can happen that 
more than one edge leads from one node to the another (they will differ in their label 
however). Irrespective of the elementary tableau a, each node in a BG can be assigned a 
depth, i.e., a unique distance from the root, and it also holds that all branching graphs 
have a unique lowest node (leaf) given by the elementary shape a itself. 

Computing the branching graph of a compound shape A with fcrows for some 

elementary shape a (where |A| = A^|cr|) requires 0(1 ^ ) '^'^ ) ^^^P^' ^^^^^ 

is the number of shapes n obtainable from A by regular removal of a. It can be 
estimated by (see Appendix B) 



Ai>ii>...>ifc>0 \ / 

The leading contribution is D"^ = O (A^'^) (because the first row Ai < |cr|A^) and thus 
we find that BG generation takes 

Cbg = O (A;l'^l N'') (14) 

time. We should point out that this is in general not polynomial, as it might appear at 
first glance. Since k is not independent of A^, for e.g. a square shape N = k'^ we indeed 
have Cbg = O (^s/N^ ' exp[y/N]^ . We will still profit from using the branching graph 
however, because even in these cases A^a grows much faster still and thus dominates the 
total complexity of computing the momenta (for more details see Appendix B). 

One can now use the efficient graph iteration described in Appendix A to traverse 
all paths through the branching graph, simultaneously building up the extended YT as 
we go. It is necessary to compute this both at once, because a modularized approach 
of extracting the paths first and then translating them one by one into extended YT 
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incurs an additional O (N) time factor coming from the fact that each path is of length 
N -1. 

The total asymptotic complexity (in both time and memory) achievable is therefore 
indeed determined purely by the number of YT on A 

CextendedYT = O (Nx) (15) 

How much is this superior to diagonalisation? The biggest differences occur for low- 
weight SU(n) representations {e.g. singlets), and for these, the number of all multiplets 
Nx grows slower than nx, the size of the corresponding total weight space. 

Take for instance N spin S — 1/2 {N even): there are ^ 1 ^1°^ — states 

but only 

( ' 

\N/2 J \ N/2-1 J N + 2 
gtot ^ singlets. 

In addition, the other methods incur the factor hw{x) because they need to repeated 
for higher weights, as described above. This factor, while trivial for SU(2), becomes 
increasingly important for larger n . 




6. Conclusion 



We have demonstrated how the extended Young tableaux method of calculating the 
eigenvalues of Cat, the generator of the cyclic subgroup "^tv C Sn, can be used not 
only for product spaces of fundamental SU(n) representations (associated with a single- 
box Young diagram), but for those of higher ones as well, if they are symmetric, i.e., 
correspond to single-row diagrams. 

Furthermore, since extended Young tableaux derive directly from the YT on a shape 
A, it is possible by exploiting the branching rule to speed up the computation oi Cn 
eigenvalues on V^®"^ to an asymptotic complexity O (A^a) in time and memory. 

BS was supported by the Landesgraduiertenforderung Baden- Wiirttemberg. 



Appendix A. Efficient graph iteration 

The branching graph (BG) is the key data-structure for implementing fast Young 
tableaux generation, but the form presented above is not yet sufficient to allow an 
efficient iteration over all paths through it. We need to both add some additional 
information to the it and use suitable data-structures to guide the iteration itself. 

To illustrate: Taking the BG as it is, we could for instance perform a depth-first 
iteration: we use a size — 1 (= max. depth) array c[.] to record to which child we 
descend to from the node fi[c]{d) at the depth d. In each step we then descend one 
level further down the graph until we reach the leaf and there, having found a new 
path from root to leaf, we add it to our result and backtrack to the closest node where 
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we can descend in a different direction. If we are only interested in the extended YT 
corresponding to the path we can built it as we descend down to the leaf, and store it 
instead of the path. 

The problem is however, that in all this we descend and backtrack step by step 
through the graph, which will take on average O (A^) steps and this brings the total 
complexity up to O {N Nx). 

To do better, we perform a pre-computation before the iteration itself adding the 
following information to each node /i: from /i we follow the leftmost path (Imp) to the 
leaf and, while we descend, push on a stack all the nodes with more than one child, 
because only these will be potential candidates to backtrack to (see Fig. lAll) . We also 
save the 'incomplete' extended YT containing all numbers N—d{fi), N—d{fi)—l, . . . , 1 
where d{fj,) is the depth of node /i. 

Precomputation of the triple K^, 1^, for all nodes of a BG for shape A built from 
elementary shapes a takes 

Cpre = O (iV) O (D^) = O (iV^+l) (A.l) 

where we used that the function counting the number of shapes /i C A obtainable 
from A by regular removal of elementary tableaux cr is bounded from above by A^'^ (see 
Appendix B). We see it is only linearly more demanding (in A^) than computing the 
basic branching graph . 

Let us now sketch an iteration process which uses this precomputed information. 
The ingredients are first the array c[.], already known from the naive iteration above 
and still needed to keep track of where we have descended to from the node yu[c]((i) lying 
at a depth d. Furthermore we introduce a stack K (the backtrack-stack) which will at 
all times contain the depths of those nodes on our current path, where we could descend 
in a different direction, i.e., those which have more than one child and have not yet 
been exhausted {c[d] < ^children of node fi[c]{d)). Flow control requires only a single 
'while'-loop which is repeated as long as K is nonempty. 

The loop performs the following steps 




{3} {3} 

Figure Al. (Colour online) To enable efficient graph iteration, we need to augment 
the basic branching graph with additional information (shown in blue): to each node 
we attach the potential backtracking-positions, i.e. all nodes with more than one child, 
lying on the leftmost path, descending from that node to the leaf. 
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0. assume we enter the loop with c[.] initiahsed and a complete extended YT E (from 
initialisation or the previous pass) 

1. first, retrieve the uppermost element (depth) from K (removing it in the process) 

2. if that element is, say, j, increment c[j] by one, reset c[i] = 1 for all i > j and 
descend to z/ := fi[c]{j + 1)) (the next, still unvisited child of node yu[c](j)) 

3. update stack K: if c[j] < ^children of /u[c](j) — )■ push j back onto K 

4. in any case: push all nodes (depths) indicated in the backtrack-list of node u onto 
K (the blue lists in Fig. |AlD 

5. obtain next extended YT: drop the indices 1, . . . , j from E and join the remainder 
with the (incomplete) YT T,^ which was added to node u during the pre-computation 

For standard YT, where the number of boxes \a\ in the elementary tableaux is equal 
to one, joining two parts of an extended YT can be done in a single step: we need only 
check wether row(j) < row(j+l). If so, we merge the leftmost column of the remainder 
with the rightmost column of T,y. If not, we simply concatenate. If we are considering 
tableaux with \a\ > 1, joining requires more steps, but can always be achieved in O {\a\) 
time. 

In all, the algorithm sketched above needs O (1) (or O (|cr|) i.g.) steps to generate 
one path/extended Young tableaux, and thus the full iteration requires O (|o"| Nx) time 
and O {N Nx \a\) memory if we store the complete list of extended tableaux. If we only 
keep the momenta, O (A^a) memory will suffice. 

Appendix B. Branching graph size 

The time required to generate the basic branching graph for a shape A built up from 
elementary tableaux a as well as augmenting it in preparation for efficient iteration 
is determined mostly by its size, i.e. the number of its nodes. This in turn is just D^, 
defined as 

D"^ = ^diagrams n with C A and fi obtained from 

A by regular removal of one or more shapes a (B.l) 

Our goal is now to find a good estimate for D^. 

Defining Dx := D'l~^^\ we can use it as an upper bound on D^, as the additional a- 
dependent constraints in (IB. II) serve only to decrease the number n that are compatible. 
But Dx is easily expressed as the multiple sum 

«A=E E ■■■ El (B-2) 

ji<Ai i2<Min(A2,ii) jfc<Min(Aft,jft_i) 

where k is the number of rows in A. This can be estimated from above by forgetting 
about Min(Ai,ji_i) and bounding ji just by instead 



Dx< 




(B.3) 
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Two instances are of particular interest (set a = 1): 

• Ai = A^, A; = 1 in which case the estimate gives almost the exact result {Dx = N 

f N+l\ 
compared to I ^ I = A* + 1) 

• Ai = N/k =: m,k > 1, A;|A^(rectangle shape) where the above estimate gives the 
exact result, as we will show in the following. 

Building the branching graph for a rectangular shape removing in turn 
1, 2, . . . , s, . . . , A^ — 1 boxes is equivalent to building (box by box) shapes with no more 
than k rows and m columns. Without restrictions, the number of shapes with N boxes 
is simply p{N), the number of integer partitions of N. With the restrictions, we must 
instead use p<k;m{n), the number of integer partitions using at most m summands 
of size < k. Thus we obtain the exact branching graph size for a rectangular shape 
A = (m, . . . , m), if we sum this over all steps s = 1, . . . ,N: 

N 

D(m,...,m) = ^P<k;mis) . (B.4) 
s=l 

However, a little thought reveals that this and the sum ( IB. 31) are, in fact, the same, 
proving that in this case the bound flB.2|) is tight and the number of nodes is exactly 



given by (m + k)\/k\m\. 

What is the asymptotic complexity in terms of A^? If <^ Ai ~ A^ (or vice 
versa), clearly (Ai + k)\/k\\i\ < {N + k)\/k\N\ = 0{N^) and therefore polynomial in 
A^. However, if k ^ y/N {e.g. shapes of square or triangular form like (/c, /c — 1, . . . , 1)), 
then our upper bound is 

and since we have shown that it is tight in the case Ai = N/k, we see that there 
are indeed shapes for which computing the branching graph is of nearly exponential 
complexity. But exactly these shapes also have the highest number of YT, growing like 
O (exp[A^]), i.e., fully exponential. Therefore computing the branching graph is always 
worthwhile, as it is in all cases much less costly than generating the YT from it. 

For many other combinations of A^ and k ( IB. 21) overestimates the size of the the 



branching graph considerably. Take A = (A^ — A; + l, !,...,!). Assuming N — k + l> k 
the true value is Da = (Ai - 1)^ + (Ai - 1) (Ai - A; + 1) ^ 2A? - A;Ai = O {X\) independent 
of k (as long as it remains smaller than Ai) while ( IB. 21) yields O (Aj). 

Appendix C. Momentum offset 6o 

Given a tableaux T built from A^ elementary tableaux a we may interpret it as pertaining 
to the product space V®^ of any SU(n) where n is at least as large as the number of 
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N = 


> Ct 


total 


half-integer 


even = 


> odd 


integer 


)) 


odd = 


> even 


77 


integer 


even = 


> even 


77 


)) 


odd = 


> odd 


57 



Table CI. The four possible parity combinations of J^c (*)c ''^^'^ ^ead to an 

integer value for the total momentum (jC.2[) . 



rows in T. We want to show here, that the momentum assigned to T via the sum ([7]) 
is independent of this interpretation, i.e., independent of n. 

J2b. = - - ^) = ^ (E <^)c - ^r) - IW\N' (C.l) 

ceE(T) c \ c J 

where we defined ct as the number of columns of E{T) and used the relations ct = J2c 
Ec^c = W\N and Ec^c {i)c = |cr|A^(A^ + l)/2. As a reminder, kc is the number of boxes 
in column c of the extended tableaux E{T) and IcjiV is just the total number of boxes 
in E{T) (and therefore also in T). 

We see, that if we add the offset momentum number bo = —{n — l)N^/2 we arrive 

at 



ceE(T) \ c 




^_W.-o(cT+|cr|iV^) I (C.2) 

and thus n cancels when computing the momentum pt = 2n/Nn{^^bc + bo). 

What still needs to be checked is whether the quantity in parenthesis in (lC.2p is 
always an integer. To see that this is indeed the case, we need to analyse the relationship 
between A^, ct and J2c (^)c- Since in all columns of E(T), the boxes are in sequence, 
Ec (*)c always either integer or half-integer. In fact we can express each summand 
as (^)c ~ + {kc — l)/2, where jc is number in the uppermost box. Therefore, (i)^ is 
half-integer, if and only if there is an even number of boxes in column c. Now assume 
Ec (^)c is half-integer. This means, we must have an odd number of columns with an 
even number of boxes. If now is even, there is an even number of boxes left to be 
distributed over rows with an odd number of boxes in them. This means that this num- 
ber of (odd-box- number) columns must be even. Thus in this case cr/2 is half-integer 
while A^^/2 is integer and in total sum flC.2p is of the form nx integer. It is not hard to 



see that in the other three cases this holds as well (see Tab. [CT 
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